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TECHNICAL  NOTE  2k32 


TRANSFORMATIONS  OF  THE  HODOGRARH  FLOW  EQUATION  AND  THE 
INTRODUCTION  OF  TWO  GENERALIZED  POTENTIAL  FUNCTIONS 
By  Luigi  Croccol 

SUMMARY 


It  has  been  shown  that  the  hodograph  equations  of  motion  can  be 
derived  in  a  symmetrical  form  by  the  choice  of  the  velocity  and  the  mass 
velocity  as  independent  variables.  The  equations  obtained  by  the  use 
of  the  velocity  potential,  the  stream  function,  or  their  transforms  as 
the  unknown  function  are  of  the  same  general  form  and  therefore  can  be 
treated  in  the  same  manner. 

Bartlcular  sets  of  solutions  have  been  studied  Independently  of  the 
gas  law  adopted  and  some  properties  of  the  series  obtained  by  means  of 
these  sets  have  been  discussed.  Approximate  gas  laws  for  which  the 
solutions  of  the  hodograph  equations  can  be  easily  found  have  been 
briefly  discussed. 

The  equations  have  been  further  transformed  so  as  to  have  as 
Independent  variables  the  complex  velocity  and  the  complex  mass  velocity. 
Two  new  generalized  potential  functions  can  then  be  Introduced  that 
satisfy  very  compact  equations.  From  these  functions,  all  the 
quantities  concerned  with  the  representation  of  the  motion  can  be 
derived  by  means  of  formulas  independent  of  the  gas  law  adopted.  By 
means  of  tke  generalized  potential  functions  some  developments  have  been 
performed  with  the  approximate  Chaplygin-Von  E^rmlin-Tsien  law. 

An  approximate  transonic  method  has  also  been  suggested. 


INTRODUCTION 


From  a  purely  mathematical  point  of  view,  the  ordinary  hodograph 
equations  for  the  stream  function  or  for  the  velocity-potential  function 
and  the  equations  relating  them  to  the  physical  coordinates  are 
sufficient  for  the  study  of  two-dimensional  Isentropic  flows.  However, 
from  a  more  physical  point  of  view  they  are  not  very  elegant  because  of 
their  lack  of  symmetry  in  contrast  with  the  symmetry  of  the  corresponding 
relations  for  the  incompressible  case. 

^At  present  at  Guggenheim  Jet  Propulsion  Center,  School  of 
Engineering,  Princeton  University,  Princeton,  New  Jersey. 
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Now^  the  eq^uations  that  define  the  relocity  potential  0  and  the 
stream  function  ^  are 


0n  =  0 

^ 


where  the  subscripts  denote  the  differential  quotients  with  respect  to 
the  element  of  streamline  ds  or  the  element  of  normal  dn,  obtained 
from  ds  by  a  counterclockwise  rotation  of  90°^  pi^)  y 

Pq  =  p(0)  represent^  respectively^  the  velocity,  density,  and  stagnation 

density.  Equations  (l),  which  are  symmetrical  with  respect  to  ds  and  dn 
if  p  =  p  ^  conserve  their  property  of  symmetry  for  variable  p  if  the 


^0^ 


mass  velocity  m  is  considered  in  some  way  the  counterpart  of  w.  If 
the  hodograph  equations  for  0  and  t  (or  for  other  functions)  can  be 
expressed  so  as  to  make  w  and  m  (instead  of  the  relation  connecting 
them)  appear  explicitly,  the  equations  will  then  have  a  symmetrical 
form  that  can  be  interesting  not  only  from  a  formal  point  of  view  but 
also  from  the  fact  that  it  can  give  rise  to  many  possible  developments, 
some  of  which  are  illustrated  in  the  present  paper.  In  particular,  it 
is  possible  to  choose  as  new  independent  variables  the  complex  velocity 
and  mass  velocity  and  to  introduce  a  new  generalized  potential  function 
satisfying  a  very  compact  equation  from  which  0  and  their  Legendre 
transforms  X  and  cn,  and  the  physical  coordinates  x  and  y  can  be 
deduced  by  simple  differentiations. 
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HODOGRAPH  EQUATIONS 


The  hodograph  equations  can  he  directly  deduced  as  follows.  If  N 
and  S  (fig.  l)  are  the  normal  and  the  suhnormal  to  the  streamline 
and  9  is  the  direction  of  motion  at  a  point  P 


z  =  X  +  ly  =  e^®(S  +  IN)  (2a) 

dz  =  ei®[dS  -  N  de  +  l(dN  +  S  d0)^ 

=  ei0(ds  +  1  dn)  =  i  W 

where  the  defining  equations  (l)  have  been  used  in  the  last  step.  It 
follows  from  equations  (2)  that 


d0  =  w(dS  —  N  d0) 

di|r  =  m(dN  +  S  d0) 


(3) 


j 

These  are  two  relations  between  exact  differentials  and  therefore  can  be 
written  as 


dS  =  dw  +  S0  d0  =  I  dw  +  0Q  +  n) 
dN  =  N^  dm  +  N0  d0  =  I  dm  +  (|  -  s) 


d0 


W 


y 


V 


k 
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vhere  m  and  w  are  two  related  varlatles  so  that  the  meaning  of 
partial  differentiation  with  respect  to  m  is 

S _ ^  ^  ^ 

5m  dm  5w  m*  5w 


Since  dS  and  dW  are  exact  differentials ^  it  follows  from 
equations  (4)  that 


lii  *  'l  .  Sjl  *  S 

50  Vm  5m\m 


so  that 


=  00 


m2Sm  =  -% 


But  from  equations  (4) 


=  i|r. 


mW-w  =  ilf. 


wSy  -  ^  or  ~  01 
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% 


Hence ^  putting  Ijj:  equal  to 

transformation  for  m  results  in 


and  performing  the  same 


the  follovlng  equations: 


(6) 


which  yield  the  well-known  Chaplygin  equations  for  0  and  t  in  a 
symmetrical  form. 

It  is  seen  from  equations  (4)  that 


and 


and,  with  the  aid  of  equations  (5), 


Sg  =  wN^  +  W  = 

bw 

and 

-^0  =  mSjjj  +  S  = 

Sm 


V 
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These  eq.uations  can  he  satisfied  hy  putting,  respectively. 


S  =  ,  wW  =  Xg  ,  w  =  c%  ,  mS  =  -CCQ  (7) 

Eq^uations  (7)  a^re  consistant  if 


(8) 


a  symmetrical  system  of  equations  in  X 
equations  (6).^ 


and  CD  very  similar  to 


From  equations  (3)  and  (7)  it  is  deduced  that 


d0  =  w  dX^  -  Xq  d0  =  d(vX^  -  X) 


and 


dilf  =  m  doim  -  ixg  d0  =  -  o)) 


Hence,  to  within  an  unessential  constant. 


0  =  wX^,  -  X 

;  (9) 


t  =  mCD^  -  CD  j 

*3 

which  give  0  and  ^  in  terms  of  X  and  cd.  The  functions  X 
and  CD  are  of  course  the  Legendre  transforms  of  0  and  ^  considered 
as  functions  of  the  physical  coordinates. 


The  functions  X  and  cd  are  distinguished  hy  the  fact  that  once  a 
solution  of  equations  (8)  is  known  all  the  other  functions  concerning 

Equations  (8)  have  already  "been  written  in  the  present  form  hy 
Bateman  (reference  l) . 

3Relations  (9)  already  have  been  derived  in  the  present  form  hy 
Bateman  and  Peres  (references  1  and  2). 
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the  physical  representation  of  motion  can  be  derived  by  simple 
differentiations^  0  and  t  being  obtained  from  equations  (9)  and  z 
from  equations  (2)  and  (7),  However ^  if  0  and  ^1;  are  the  known 
functions  satisfying  equations  (6)^  integrations  are  necessary  to  deduce 
the  other  quantities.  After  integration  and  determination  of  the 
constants  so  as  to  satisfy  equation  (8)_,  the  following  explicit 
expressions  for  X  and  o)  are  obtained  from  equations  (9): 


X  = 


0(wr,0l)3ln(e  - 


1 


t  (wp^^l)cos(^0  —  02)d0^  +  cos  0  +  C2  sin  ^ 


03  =  -HU 


^Cwq^ejd 


0(w^^0jl)co3(0 


0l)d0i  —  sin  6  +  cos  0 


where  w^  and  0^  are  two  arbitrary  reference  quantities,,  m^  =  m(wr)^ 
and  and  C2  are  two  arbitrary  constants  with  no  influence  on  0 

and  It  is  readily  deduced  with  the  aid  of  equations  (2)  and  (7) 

that 


-  i 


i0T  , 

e  -^-de^  -  (Cl  +  IC2) 


By  differentiation^ 

V 


which  agrees  with  equation  (2b). 
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From  equations  (6)  and  (8)  the  hodograph  equations  for  the 
functions  x>  to,  0,  and  ijr  can  he  easily  shown  to  be  as  follows 


w  —  m  —  I  +  =  0 

8m  \  dw/ 


+  ^=0 


(10a) 


(10b) 


l8 _ Pi  80  1  8^0  ^ 

“Ki) 


l8 _ fl  dir]  8^t  ^ 

^h(l\  ^802" 


(10c) 


(lOd) 


Each  of  these  equations  reduces  to  the  Laplace  equation  if  m  =  w. 

The  following  equations  are  obtained  from  equations  (lO)  with  w 
as  independent  variable: 


v2  ^  t  n  +  m2  _  V  i-  log  (l  -  M2)1w  1“  +  (l  -  ^  .  0 

7^2  I  dw  ®  8w  ^  ^  p.fl2 


w2  ^  +  n  _  m2  _  w  —  log  (l  -  M^j]  w  M  +  (l  -  ^  =  0 

8w2  L  dw  «  '  8w  ^  "  .,2 


(l  +  M2)w|i+  (i_m2)^=0 


2 


NACA  TN  2^32 


9 


where  M,  the  Mach  nimiber,  is  defined  hy  (see  equation  (15)) 


V 

1  _  ^  ”• 
d  log  w 


For  every  particular  law  m(w);  that  is,  for  every  p(w),  explicit 
equations  are  obtained.  For  M  <1  (subsonic  flow)  the  equations  are  of 
the  elliptic  type;  for  M  >  1  (supersonic  flow)  the  equations  are  of  the 
hyperbolic  type. 

The  holograph  equations  are  frequently  transformed  so  as  to 

simplify  the  second-order  terms.  Thus,  if  for  -3—  >  0, 

CLv 


then 


^  =  a^dX 

V 


dv 

m 


m  d(i)  =-—  dX, 

wy  32 


V 
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For  —  <  0  analogous  transformations  can  te  performed  vith  the 
dw 

Introduction  of  dX^^,  a^,  and  defined  in  the  same  way  as  dX,  a, 

and  P  with  -dm  instead  of  dm.  It  is  seen  that  X  represents  some 
kind  of  an  integral  mean  between  log  w  and  log  m.  Introducing  the 
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1/2 


and  the  Mach  numher  M  =  — 

a 


souQd  velocity  a‘= 
follovin^  equations  deduced  from  the  Bernoulli  equation 


leads  to  the 
dp  +  pw  dv  =  0: 


\ 

a2  +  v2  =  0 

d  log  p 


il2§-£  =  ^2  ) 

d  log  w 


d  log  m  P 

- =  1  _ 

d  log  w  / 


(15) 


Then  for  M  <  1, 


dX  = 

\/l  -  d  log  w 

(16a) 

“H 

^PO  ) 

(l6b) 

(l6c) 

There  is  obtained  a  corresponding  set  of  equations  for 

dXj^  j  ^  and 

when  M  >  1  by  simply  replacing  1  —  by  —  1. 

Approximate  Methods 

It  is  seen  from  equations  (I3)  that  for  constant  a  or  p  these 
equations  reduce  to  the  Laplace  equation.  The  first  possihilitj  is  to 
he  rejected  for  subsonic  motion  because,  as  equation  (16)  shows,  it 


12 


MCA  TN  2432 


gives  p  as  an  increasing  function  of  M.  The  second  one  is  the  well- 

known  Tsien  approximation  (references  3  and  4)  t 


=  —  n/l  -  =  'K 

P 


(17) 


where  K  is  a  constant.  This  equation  reduces  to  the  Chaplygin 
approximation  for  K  =  1  (reference  5  and  derived  studies). 

Some  considerations  that  will  be  useful  in  a  subsequent  section  are 
now  introduced.  With  the  help  of  equations  (12)^  equation  (17)  can  be 
immediately  integrated  to  obtain 


where  H  is  an  arbitrary  constant.  Integrating  now  equation  (l6a)  and 
the  Bernoulli  equations  (15)>  there  follows^  respect ively^ 


X  =  log 


w 


+  \/l  +  Hw^ 


+  Consteint  =  — logfi  +  -i- ]  +  Constant  (199') 
V  rnSfKj 


and 


0 


1  pq 

-  +  Constant 

EH  p 


(19h) 


MCA  TN  2432 


13 


Prom  equations  (I5),  for  a  real  Isentropic  gas. 


^  =  /l  -  ^ 

-2  „  2 


1  +  -21 - - 


(20a) 


a^ 


(20b) 


where  w^^^^  _  aQ^ - ^  Is  the  maximum  velocity  at  zero  density, 

a^  is  the  stagnation  sound  velocity,  and  7  la  the  adiabatic  index. 

A  comparison  between  equations  (18)  and  (20)  is  shown  in  figure  2 
where  the  law,  equation  (I8),  is  represented  by  a  straight  line. 
Chaplygin  takes  it  as  the  tangent  to  the  graph  of  equation  (20) 

at  w  =  0,  so  that  K  =  1  and  H  =  -i-.  In  the  Karm^-Tsien  method^ 

ao2 

the  tangent  is  taken  at  w  =  w^^,  the  speed  at  infinity,  so  that 


KE  = 


2 

PQ  1 

2  2 

P«. 


a  (1  -  M  ^ 
00  \  00 


^This  method  has  been  often  presented  in  a  less  coherent  form,  as 
the  constants  of  equations  (l8)  and  (I9)  are  determined  for  different 
conditions^  though  the  formula  for  the  correction  of  pressure  coefficients 
is  not  affected  hy  this  incoherence. 


li^ 
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Observe  that  for  K  1  the  value  of  p  at  w  =  0  differs  from  the 
exact  stagnation  density  Pq.  It  is  easily  seen  that  the  value  for  KH 

in  equations  (21)  gives  the  slope  of  the  true  isentropic  relation  at  a 
point  corresponding  to  the  conditions  at  infinity j  hut ,  as  this  depends 
only  on  the  value  of  KH,  it  is  seen  also  that  K^Lrman’s  condition  is 
satisfied  for  every  parallel  to  the  tangent  at  v  =  Va,  (fig.  2),  that 

is,  for  the  KH  value  given  hy  equation  (21)  hut  for  different  values 
of  K.  This  suggests  the  posslhility  of  improving  the  mean  approximation 
of  the  Kann^— Tsien  method  hy  an  appropriate  choice  of  K.  The  KannAn— 
Tsien  correction  formula  for  the  pressure  coefficient  has  then  to  he 
modified.  The  modified  formula,  independent  of  the  constant  of 
integration  in  equation  (l9a),  is  then 


which  reduces  to  the  Karmm-Tsien  correction  formula  for  the  value  of  K 
given  hy  equation  (21).  Application  of  the  modified  formula  with  some 

2 

value  of  K  between  —  (l  -  )  and  1  gives  values  of  Cp  in 

2 

Poo 

better  agreement  with  experimental  values. 

For  supersonic  motion  the  hodograph  equations  reduce  to  the 
simplest  hyperbolic  equation  (wave  equation)  for  =  Constant 
oj.  =  Constant with  and  given  by  equations  (12)  with  -dm 

Instead  of  dm  or  hy  equation  (l6)  with  —  1  in  place  of  1  —  M^. 
The  second  possibility  is  now  to  he  rejected  because  the  resulting  value 
of  p  increases  with  M.  The  first  possibility  gives 
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that  iSj  after  integration 


5 


w2  = 


The  constants  and  can  he  determined  so  as  to  satisfy  E6nni^n*s 

condition; 


%% . 


a 


4 

GO 


H. 


aoo%j^ 

Moo^  -  1 


With  =  Constant j  equations  (I3c)  and  (l3d)  (modified  for  M  >  1) 

reduce  to  the  simple  vave  equation.  The  general  solution  can  therefore 
he  represented  hy,  for  Instance, 


X  =  r(xi  ±  e) 


(22) 


5 Peres  has  already  indicated  a  law  of  this  kind  (reference  2). 
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Other  laws  for  which  the  holograph  equations  reduce  to  the  Laplace 
equation  for  M  <  1  (or  to  the  wave  equations  for  M  >  l)  are  imme¬ 
diately  deduced  from  equations  (li<-)^  (or  from  the  corresponding  equation 
for  the  supersonic  case)  as  the  laws  which  make  one  of  the  four  quan¬ 
tities  a,  l/a,  P,  and  l/P,  (or  those  corresponding  for  M  >  l)  linear 
in  X.  (or  in  X]_).  For  the  true  Isentropic  gas  the  curves  of  a,  P, 
l/a,  and  l/P  as  functions  of  X  are  shown  in  figure  3.  For  w  =  0, 

X  =  —CO,  and  a  =  P  =  1  the  Chaplygin  and  the  KarmAn-Tsien  approximations 
replace  the  true  shape  hy  a  horizontal  line  which  can  give  an  approx¬ 
imation  not  too  had  even  for  w  =  0.  Approximate  laws  for  which 
equations  (l4)  reduce  to  the  Laplace  equation  are  represented  by 
arbitrary  (and  generally  not  horizontal)  sti^i^t  lines.  It  is  seen 
that  for  X  =  —CO  these  lines  diverge  hopelessly  from  the  true  law. 

Hence  these  laws  do  not  appear  to  be  convenient  for  the  approximate 
representation  in  a  large  range  of  velocity.  Nevertheless  they  can 
possibly  have  application  when  the  variations  of  velocity  from  a  mean 
value  (for  Instance,  the  value  at  Infinity)  are  small.  In  this  case  it 
is  possible  to  achieve  a  better  approximation  than  with  the  EarmAn- 
Tsien  method  by  taking  as  the  approximate  law  the  tangent  to  one  of  the 
curves  of  figure  3,  The  resulting  approximate  p,p  curve  will  have  a 
contact  of  second  order  with  the  real  isentropic. 

For  the  supersonic  case,  a]_,  P2_,  l/n-p,  and  I/P^j  fi^nctlons 
of  X^,  are  shown  in  figure  4  for  7  =  1.4.  An  interesting  possibility 
is  given  by  the  curve  of  a^,  which  can  be  well  approximated  by  a 

stral^t  line  between  M  =•  2  and  M  =  10.  Hence  in  this  range  the 
exact  hodograph  equation  in  =  a  X  differs  very  little  from  the 

simple  wave  equation  in  X^  and  0.  Therefore  the  general  solution  of 

the  supersonic  motion  in  the  said  range  of  M  is  approximately 


X  =  ^  f(Xi  +  9) 

tti 


where  and  are  given  hy  the  true  isentropic  law.  This 

approximation  seems  to  he  better  than  the  approximation  given  hy 
eq^uatlon  (22). 
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Exact  Solutions  of  The  Hodograph  Eq^uations 

The  pQ-yrer  set.—  Many  studies  have  been  developed  on  two  sets  of 
particular  solutions  of  the  hodograph  equation  in  0  or  the  so- 
called  power  set  and  exponential  set,  characterized  by  the  fact  that 
for  the  incompressible  case  they  reduce  respectively  to  the  natural 
powers  of  the  logarithms  of  the  complex  velocity  and  those  of  the 
complex  velocity  itself.  The  symmetrical  form  of  eqxiatlons  (lO)  makes 
it  possible  to  present  these  solutions  in  a  form  that  seems  interesting. 

By  the  introduction  (thou^  it  is  not  strictly  necessary)  of  four 
auxiliary  quantities  cS,  %  and  satisfying  in  the  same  order 
equations  (lO),  and  some  supplementary  conditions,  the  four  complex 
quantities 


F  =  X+  i(D 


F  =  X  +  Ice 


G  =  0  +  it 

and 

S  =  0  +  it 


can  be  made  to  satisfy  not  only  equations  (10)  in  the  same  order  but 
also  the  relations 


mF  = 

(23a) 

li 

(23t) 

1 

V  ^1/m  - 

:  -iGg 

(23c) 

m  " 

-IGq 

(23cL) 

and 
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corresponding  to  equations  (6),  (8),  and  the  relations 


G  =  wF,  -  F 

W 


jf\f 

a  =  jnF^  ^  F 


{2k) 


corresponding  to  eq^mtlons  (9).  Conyersely  F  and  F  are  given  by 
formulas  corresponding  to  those  -written  in  section  entitled  "Hodograph 
Eq^uations”  for  x  and  o): 


F  = 


G(w,0)d(± 


Vn 


G(w^,ei)sin(0  -  0i)de3^ 


+ 


G(v2;,j0i)cos(0  —  0]_)d0i 


(25) 


and  a  similar  eq^uation  for  F  vith  m  and  -w-,  G  and  G  interchanged. 

The  quantities  X  ^  ^  determined  so  that  they  reduce 

to  X,  G^y  ^y  and  for  the  incompressible  case  m  =  -w;  hence  the 
equations 


Fi  =  Fi  =  Xj^  +  ±(D^ 

Gi  =  Gj  =  +  Itj 


and 
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satisfy  the  Laplace  equations 


AF^  =  0 


=  0 


where  A  -  -  +  - |  to  which  equations  (lO)  reduce  when  m  =  w, 

V  s(iog  w)^ 


An  operator  ((  ))  is  now  defined  as  that  which,  when  applied  to 

.h 

log  — )  (w^  being  an  arbitrary  reference  velocity),  transforms  it 


dwi  (im2  p^2  dW3 


2 


(26a) 


for  h  =  1,  2,  .  .  the  integration  heing  repeated  h  times;  then,  the 
indication  of  the  lower  limit  of  integration  being  omitted  for  brevity, ° 


^r\h 


“  dm-i  p^l  dwo  ^“2  dm 


m 

r 


W^  J  22]^  ^3 


^  =  h'l(m,w)  (26b) 


mi  d| 


ji/m^  •  •  •  =  (26c) 


m-r.\h 


l/nv  wy  i/„ 


formulas  (26c)  and  (26d)  could  preferably  be  written  as  operations 
iwern  of*  1  ogt  iilr/y  and  log  -t _ 


on  powers  of  log 


and  are  different  from  the 


operations  given  by  equations  (26a)  and  (26b)  on  powers  of  -log  ^ 

^r 

and  -log 

iiip 
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From  these  eq^uations  there  is  deduced 


and  the  analogous  equations  obtained  by  interchanging  v  and  m_,  and/or 
by  changing  w  into  l/w  and  m  into  l/m. 

It  is  then  easily  verified  that^  vith  =  1,  each  of  the  following 
functions 


F 


F 


G 


G 


is  a  solution  of  the  corresponding  equation^^  the  signs  having  been 
selected  so  as  to  satisfy  also  equations  (23).  The  choice  of  the  upper 
or  lower  sign  does  not  affect  the  values  of  0^  and  o)  derived 

from  equations  (28)  so  that  for  the  solution  of  the  flow  problem  it  is 


7  . 

‘This  kind  of  solution  has  been  first  discovered  by  Bergman  (see 
for  instance  reference  6)  and  by  Bers  and  Gelbart  (reference  7).  The 
present  form  is  new  and  more  symmetrical. 
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sufficient  to  retain,  for  instance,  only  the  upper  one;8  nevertheless  in 
some  cases  it  can  he  useful  to  consider  solutions  with  both  signs. 


The  solution  represented  by  the  values,  equations  (28),  of  G 
and  6f  does  not  coincide  with  the  one  corresponding  to  the  values, 
equations  (28),  of  F  and  F  for  it  differs  from  the  values  of  G 
and  Sf  derived  from  F  and  by  means  of  equations  (24) .  This  can 
be  shown  as  follows :  From  the  recurrence  formula,  which  la  easy  to 
verify. 


and  from  the  values,  directly  deduced  with  Iq  =  1, 


and 


l2(w,m)  -  ^  l3^(m,w) 


There  is  deduced  with  the  aid  of  equations  (27) 
_  V  -|^jljj(w,m)  =  Ij^(w,m)  -  g  Ij^_3^(m,w) 


(29) 


and  the  analogous  relations  obtained  after  interchanging  the  variables, 

Q  — -- 

In  this  case  for  m  =  w  the  four  solutions,  equations  (28),  reduce 

(to  within  a  multiplicative  constant)  to  ^log  ^ - lej^,  that  is  to  the 

power  set  for  the  incompressible  case. 


2.2 
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Application,  of  equations  (24)  to  equations  (28)  yields  the 
following  ec^uation: 


and  a  similar  equation  for  -Gj  if  w  is  Interchanged  with  m,  with 
the  ±  sign.  This  G  differs  from  the  elementarj  solution, 
equations  (28),  although  it  is  a  linear  combination  of  such  elementary 
solutions,  a  fact  holding  also  for  the  incompressible  case.  In  that 
case,  however,  the  expression  for  G  contains  only  an  elementary 
solution  with  the  upper  sign  if  F  is  so,  for  then  w^  =  m^. 

The  physical  coordinates  are  easily  deduced;  for  by  equations  (7) 


(30) 


X  and  CD 


z  =  ei®(S  +  IN) 


with 
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It  is  now  again  verified  that,  as  already  observed,  the  choice  of 
the  sign  does  not  affect  the  results  concerning  the  values  of  0,  ilr, 
and  z  hut  only  the  Introduced  auxiliary  factions .  In  subsequent 
work,  therefore,  only  the  upper  sign  la  retained. 

Linear  combinations,  and  in  some  cases  infinite  series  of  the 
elementary  solutions,  equations  (28),  are  still  solutions  of  the 
corresponding  equations. 

Infinite  series  in  the  power  set.—  If,  say,  F  is  given  by  an 
infinite  aeries,  then  developing  and  inverting  the  order  of  summations 
gives 


-  2:  £ 


n=h  ( 


T  “n(-i®) 


Ij^(w,m) 


5= -10 


■where 

0 


represents  the  function  corresponding  to  this  power  series. 

Now,  independently  of  the  convergence  of  this  power  series  it  la 
readily  verified  by  differentiation,  that  equation  (3I),  whenever  it 
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¥ 


The  formula  for  F  can  he  written  in  a  different  form: 


F  = 


(w,m) 


h=0 


^^\e)  +  iap^^^e) 


>] 


with 


A(-i0)  =  ajL(e)  +  ia2(0) 


and 


B(i0)  =  h^(0)  +  ihgCe) 


It  follows  that 


with 


and 


+  ih2^^^0)J 


h3^(0)  =  ^^(0)  _  ^a2*(0) 


^2^®)  ~  +  —  an  *(0) 

Wv» 


There  are  similar  eq.uatlons  for  F  and  G-. 
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Hence  for  w  = 


so  that  A  and  all  the  solutions  of  eq.uations  (31)  and  (32)  are 
determined  hy  the  values  of  X  and  Its  radial  derivative  on  the 
circle  w  =  Wj.  of  the  holograph  plane  or  hy  the  corresponding  values 

for  CD.  Similar  statements  hold  for  B,  and  and  solutions  of 
actuations  (33) •  Hence  the  solutions  written  depend  upon  two  arhltrary 
functions  and.  In  their  region  of  convergence,  represent  the  general 
solution  of  the  holograph  eq.uatlons.  Naturally  they  do  not  give  any 
Indication  of  the  Behaviour  of  the  corresponding  solutions  at  w  =  0; 

for  as  m-^  0,  “  as  |log  so  that  the  origin  Is 

certainly  outside  the  region  of  convergence. 

In  fact.  If  the  solutions  must  he  regular  at  w  =  0,  only  one  of 
"tHe  functions  a(0)  and  h(0)  can  he  chosen  arbitrarily,  and  other 
representations  of  the  solutions  are  needed  to  determine  the  other.  The 
region  of  convergence  of  the  series,  equations  (31),  (32),  and  (33) ^ 
depends  on  the  form  of  A(0  and  B(5).  However,  a  general  Idea  of 
Its  shape  can  he  given  hy  mahlng  very  general  assumptions  about  these 
functions  and  m(w). 

Let  r(0)  he  less  than  the  distance  In  the  plane  between  the 
point  —10  and  the  nearest  singularity  of  A(^  )  and  let  Ajj^(0) 
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"be  the  upper  hound  of  |a(^  )|  on  the  circle  with  center  —10  and 
radius  r.  Then  the  Cauchy’s  Inequality  gives 


4?'*  c.-ie  rCe)** 


(35) 


and  a  similar  expression  for  with  and  r^^  in  place  of 

A^av  and  r.  For  the  If 


(l^ _ 'l 

Vm  d  log  w^jj^  V  P  /inax 


1  dm 
w  d  log 


m 


d-(l/w) 
d  log  U/w) 


max 


>  (36) 


w 


d(l/m) 

^  log  (1/w) 


max 


are  the  upper  hounds  of  the  written  quantities  between  w  and  w^^ 
equations  (26)  give  for  even  values  of  h 


and  the  same  limitations  for  Ij^(m^w) 


and 


Ih(l/m,l/w)  . 
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For  odd  values  of  h  the  upper  hounds  for  |l2^(wjm)| 
and  are  obtained  from  the  corresponding  expressions  in 

1  /2  1  /2 

equations  (37)  hy  multiplying  them  with  (a/h)  ^  and  (h/a)  ^  , 
respectively^  and  those  for  |lj^(l/w_,l/m)j  and  llj^Cl/m^l/w)!  are 

obtained  by  multiplying  the  expressions  in  equations  (37)  hy  (c/d)^/^ 
1/2 

and  (d/c)  '  ^  respectively. 

Consideration  of  the  series  (31)  shows  that  the  terms  of  the  series 


are  less  than  the  corresponding  terms  of  the  geometric  series 


which  converges  when  the  ratio  is  less  than  1. 

Hence  the  series  (3I)  and  (as  it  may  be  deduced  in  the  same  way) 
(32)  converge  absolutely  in  the  region 

(at)V2  |iog  <  r(0)  (38a) 

Similarly^  the  series  (33)  converge  absolutely  for 


(cd)^/^  jlog  <  rj^(a)  (38-5) 

Nov  since  at  and  cd  are  functions  of  w  and  and  since  r 

and  are  quantities  vhlcli  increase  with  the  distance  from  the 

singularities  of  A  and  B,  the  general  shape  of  the  region  of 
convergence  in  the  holograph  plane  is  a  curved  strip,  which  contains  the 
circle  w  =  w  and  whose  width  will  he  a  minimum  when  -10  is  nearest 
to  a  singularity  of  A(^  )  or  B(0' 
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If  a  singularity  lies  on  the  imaginary  aocis,  that  is,  if  A(-i9)  or 
B(i0),  and  their  derivates  have  a  singularity  for  some  value  of  0,  the 
corresponding  width  of  the  region  of  convergence  will  he  zero.  This 
happens,  for  instance,  when  the  reference  velocity  w^,  is  the  velocity 

at  infinity  of  the  flow  round  a  body. 


An  observation  of  some  interest  is  that  as,  for  plausible  laws, 
p/pq  0(l)  for  M  =  1  then  (ab)^/^  and  (cd)^/^  (equations  (36)) 

(|i-m2|V2) 


are  0 
and  near  1 


or  0^(l  - 


(the  larger  of  the  two)  for  M 


Hence  equations  (38)  show  that  for  given  A(n  and  3(0  the 
width  of  the  region  of  convergence  is  the  greatest  and  the  rapidity  of 
convergence  the  best  near  the  sonic  line  M  =  1.  Therefore  it  is 
believed  that  the  solution  represented  by  series  (31),  (32),  (33) 

may  have  applications  in  the  solution  of  transonic  problems,  naturally 
in  combination  with  other  methods  converging  in  the  rest  of  the  field 
of  motion. 


Finally  observe  that  the  development  of  equation  (3I)  can  be 
handled  differently  so  as  to  obtain  a  power  series  in  0 ; 

h=0  h!  n=h 


which  by  use  of  equation  (27)  and  similar  expressions  and  with 

f^(w)  =  n!a^I^(w,m)  and  f2(w)  =  _ n!a^I  (m,w)  becomes 

n=0  n=0 


flQ  (2k) I 


i-ie) 


2k+l 


So  (2k  +  1) 


w 


(im 


Similarly^ 


F 


(-18)  2k 
(2k)! 


(-±9)2k+l  ^  ^ 

(2k  +  1)1  liw 


The  functions  F  and  F,  whenever  the  series  converge,  satisfy  the 
corresponding  equations  and  the  relations  (23),  as  can  be  directly 
verified,  with  arbitrary  fj  and  fp.  Analogous  solutions  hold  for  G 

and  G  with  10  In  place  of  —10,  l/m  and  l/w  in  place  of  w  and  m, 
and  two  arbitrary  functions  g2^(w)  and  gQ(w) .  For  real  values 
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of  fj_  and  f2  the  real  and  imaginary  parts  of  F  and  ^  are 
olDtalned  directly;  so  that  the  equation 


depends  only  on  correspondingly,  and  g^  hecome  real  and 

the  equation 


depends  only  on  g^(w) . 

The  physical  coordinates  are  then  found  hy  means  of  equations  (2) 
and  (7)  to  he 


MCA  TN  2432 


31 


From  X  +  4a  there  follows  hy  means  of  equation  (9)  an  expression 
for  0  +  i\]/-  which  must  coincide  with  the  written  one  if 


«!<">  ■  ("  L  - 


or 


From  this  coincidence  the  following  interesting  foimiulas  can  he  deduced: 


which  are  easy  to  verify  directly. 


The  meaning  of  the  written  solution  is  readily  found  hy  observing 
that  when  the  series  converge,  for  9=0,  then  x  =  f^Cw),  0  =  gj^(w), 

n  n 

(0  =  y  =  0,  X  =  — ,  and  y  =  0,  so  that  f^  =/x  dw,  g^  =  dx,  and  the, 


whole  solution  is  determined  when  the  "axial”  law  of  distribution  of 
velocity  is  given. 


Hence  the  solution,  under  a  somewhat  different  and  more  explicit 
form,  reduces  to  the  one  studied  by  Llghthlll  (reference  8)  in  his 
work  on  the  transonic  flow  in  symmetrical  channels.  As  Lighthill 
observed  the  coefficients  of  the  expansions  become  infinite  at  sonic 
speed  (for  then  dm  =  O),  so  that  the  aeries  diverge  in  the  transonic 
region.  In  this  region,  however,  the  solution  can  be  found  by  following 
the  Llghthlll' 8  ingenious  method,  that  is.  Inverting  the  series 
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clvlnG  ilk  or  the  one  giving  - - — - — >  since  the  coefficients  of  the 

(i0) 

inverted  series  are  finite  at  sonic  speed.  The  application  of  Lighthill's 
method  can  he  made  easier  hy  the  present  form  of  the  solution. 

The  exponential  set. —  If  in  e(iuatlons  (Sl)^  (3^)^  and  (33) 


A(0  =  = 

nii 


and 


B(0  = 


n 


with  arbitrary  n,  then  the  fimctions 


F  =  Ej^(w,m)e 


-InS 


F  =  E^(m,w)e 

nVm-’w/ 

« -  ^Jk’B 


— in9 


Ina 


ine 


with  coefficients  defined  hy 


E  (w,m)  =  ^  n%j^(w,m) 
h=0 


uu 

En(m,w)  =  Y1  nlilji(m,w) 
h=0 


(39) 


m 
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and  similar  equations  for  the  other  'be  solutions  of  the 

corresponding  equations  (lO)  and  (23)  vhlch  reduce,  for  m  =  v 

to  (we~^^/vj^  (equations  (lOa),  (lOb),  (23a),  and  (23b))  and 

to  (equations  (lOc),  (lOd),  (23c),  and  (23d)), 

that  is,  to  the  exponential  set  for  the  incompressible  case. 

It  is  Immediately  verified  by  means  of  equation  (37)  that  the 
series  converge  for  all  values  of  v  and  for  which  a,  b,  c, 

and  d  are  limited,  that  is,  for  which  p  and  M  are  limited  and  not 
zero.  For  plausible  gas  laws  this  excludes  only  the  values  w  =  0 
and  w  =  w^„_. 


It  is  therefore  seen,  and  easily  verified  directly,  that 
equations  (40)  and  the  two  other  equations  for  E^  are,  in  order, 

solutions  of  the  ordinary  differential  equations 


w  — im  — 
dm\  dw 


-  -  =  0 


(4la) 


d  /  <4% 

m  —  w  - — 
dw  \  dm 


-  n^<u^  =  0 
n 


(4lb) 


(4lc) 


V  d(l/m)[m  cL(l/w) 


-  n  =  0 


(4ld) 


deduced  from  equations  (lO)  by  talcing  x,  o),  0,  and  i  as  the  product 
of  a  sinusoidal  factor  in  nO  and  of  the  corresponding  function 

0n-»  normal  isentroplc  law,  equations  (4l) 

become  the  known  equations  of  the  hypergeometric  type  and  have  been  the 
object  of  the  Investigations  of  many  authors.  (See,  for  instance, 
references  5,  6,  8,  9,  10,  11,  and  12.)  Equation  (W)  has  generally  b; 
studied  with  particular  regard  to  those  solutions  that  satisfy  the 

condition  that  unity  at  w  =  0. 
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It  is  immediately  verified  that  the  solutions  do  not  satisfy 

this  condition  and  that  the  corresponding  conditions  at  v  =  Wp  are: 

(En)^  =  1 


dw 


'“•Jr  ■  ["  S  ®n(">“)]r  '  [l  dlW 


m  d(l/v)  ^n(v"iii)  r  ^ 


The  symmetrical  form  of  eq.uatlons  (4l)  allovs  some  general 
relations  to  he  easily  derived.  Some  of  these  relations^  that  may  he 
useful  for  further  developments^  are  now  stated  briefly. 


The  can  he  considered  as  the  superposition  of  two  Independent 

_ 9  — p 


solutions^  of  equations  (4l): 


Cn(v,m)  ='^  n2ki2^(w,iri) 
k=0 


S^(w,m)  =  ^  n^^+ll2^^3_(w,m) 
k=0 
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for  which  the  conditions  at  w  =  Wv.  are: 


-)  = 


f&C. 

l-dF'  =0 


(Sn)r.  =  0 


^  ^  /  dSn(m,w)\  ^  ^ 

'r  V  dm  y-r 


dw 


They  are  connected  hy  relations  similar  to  those  connecting  the 
exponentials  and  the  hyperholic  cosine  and  sine: 


C_^  =  C 

n 


'  -«n 


2C„  =.  %  +  E 


-n 


-  K„  - 
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Cn(w,iii)Cn(iii,v)  -  Sn(w,m)Sn(m,w)  =  i 
En(w,m)E_n(m,w)  +  En(m,-w)E^(w,m)  =  2 


(iCn(v,m)  , 

m - - - -  nSn(m,v) 

dw 


dSn(v,m) 

m  - — - nCn(m,v) 


m 


dv 


Similar  relations  hold  vhen  w  and  m  are  replaced  by  their  inyerse 
values. 

All  the  solutions  of  equations  (4l).  can  be  represented  by  linear 
combinations  of  and  but  these  functions  (as  all  the  series 

in  l^i)  are  not  suited  to  give  the  behaviour  of  solutions  near  v  =  0. 
Since  this  behaviour  is  very  important  for  many  physical  applications^ 
it  is  necessary  to  follow  a  different  method  of  investigation:  Let 


F  -  (Fne”“)" 

G  - 

G  = 


and 
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be  solutions  of  the  corresponding  equations  (lO)  which  reduce  to  the 
exponential  set  for  the  Incompressible  case.  Then  X  Y 

aJid  are  solutions  of  equations  (4l) .  Now,  If  for  negative  values 

of  n  the  last  two  must  coincide  with  0^  and  w  (such  that,  for 

Instance,  —  =  1  at  w  =  0)  and  If  for  positive  values  of  n  the 
first  two  must  coincide  with  the  corresponding  and  it  Is  seen 
that  Xjv,  Xj^/m,  riY^,  and  must  be  equal  to  unity  at  w  =  0.  It 


Is  now  shown  that  this  Is  possible  for  all  values  of  n 
exceptional  values:  Let 


except  some 
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and  must  satisfy  the  corresponding  eq^uatlons 


and  that 

Xi-'  XX 

with  m  and  w  interchanged  or  Inverted.  Furthermore^  these 
q,uantities  are  connected  by  the  following  relations i 


~  1 
’^n 


T  =  -= 
n  T 


(^5a) 


(45t) 


n 


nm  —  Rj^w 
”  nwkj^  —  m 


(45c) 


It  is  deduced  from  equations  (43)  that 


(46a) 


m 


(46h) 


(46c) 


(46d) 


so  that  these  quantities  are  equal  to  unity  and  analytic  at  v  =  0  if 
all  the  Integrands  are  analytic  there.  It  can  he  shown  in  fact  that. 


if 


m  _  P 

w  "  PQ 


is  an  analytical  function  of 


w  near 


w  =  0^  the  integrands 


of  eqLuatlons  (46a)  and  (46b)  are  zero  and  analytic  at  w  =  0  for  all 
values  of  n  except  negative  integral  and  half— integral-  values  (for 
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OIJJ  the  negative  Integers  If  ^  I3  „  anajljtlo  funetion  of  ^), 

while  the  Ihtegrmde  of  eiuatlone  (We)  and  (46d)  are  aero  and  analytic 
ln+.„  1°  f  values  of  n  except  positive  Integral  and  half^ 

integral  values  greater  than  1  (for  only  the  positive  Integers  ^Lur 

than  1,  if  —  la  analytic  in  w^). 


These  resets  reduce  to  the  well-knovn  results  when  the  equations 
^e  hyper^^etrle.  In  this  case  llghthlll  (reference  8)  has  g“ei  ?h, 
moot  ccmplete  aiscueslon  of  the  solutions  of  the  equation  In  md 

deduced  important  theorems,  some  of  which  may  possihly  he  Generil7er^ 

po“V-^  e^s^forL 

pole  at  n  -  -1  for  the  e<iuataons  in  0^  and  ig  a  general 

property,  which  does  not  hold  for  the  equations  in  and  o^. 

eaiiatioL'»^  solutions  of  equation  (44)  (and  of  the  analogous 

equations)  and  the  corresponding  solutions  of  equations  (24)  are:^ 


E,  =  a  =  T 

1  w  1 

(^7a) 

En  =  =  T 

1  m  1 

(47b) 

F  =  X^e-10  =  we-i® 

(47c) 

F  =  X^e-iS  =  me-i^ 

(47d) 

G  =  Yqei®  =  i  ei0 

(47e) 

G  =  Yne^®  =  i 

-L  W 

(47f) 

4-0 
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It  is  seen  from  eq.uations  (24)  tliat  the  values  of  G  and  G  (hence 
of  0  end  i)  corresponding  to  equations  (47a)  and  (47h)  are 
identically  zero.  In  fact^  it  is  seen  from  equations  (2)  and  (3)  that 
"this  solution  represents  merely  a  displacement  of  the  origin  of  the 
physical  coordinates.  Equations  (47e)  and  (47f)  are  more  interesting 
as  'they  coincide  "with  hh©  ■well— knovn  Eingleh  solu*tlon  (reference  13 )  • 
The  corresponding  values  of  F  and  F  are  determined  hy  equation  (25) 
(and  the  analogous  eciuatlons) .  ThuSj 


plus  a  constant  multiple  of  w©""^^  and  me  The  coefficients 

of  e^^  in  these  formulas^  together  with  w  and  m  respectively 
represent  two  independent  solutions  of  ©Q.uations  (4la)  and  (4lh) 

for  n^  =  1.  The  operators  — 1  +  w  ^  and  —1  +  m  ^  which  generally 

allow  the  deduction  of  two  independent  solutions  of  the  equations 
in  and  from  two  independent  solutions  of  the  equations 

in  Xn  and  suffer  an  exception  for  n  =  1  as  they  produce  only 

one  set  of  solutions;  that  is,  l/m  and  l/w,  respectively;  for  when 
applied  to  w  and  m  the  result  is  zero.  This  exceptional  case  is 
explained  in  the  section  entitled  ”A  new  set."  When  n  =  1 
and  Et  =  —  equation  (45c)  hecomes  indeterminate,  hut  the  corresponding 

value  of  T_^  can  he  deduced  from  this  relation  as  the  limiting  value 

for  n  ->  1  of  the  indeterminate  expression.  It  follows  from 
equation  (44)  that  with  the  condition  =  1  at  w  =  0 
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Similar  expressions  hold  for  and  Y_3_.  Hence  the  equations 


represent  the  second  solution  of  equations  (4lc)  and  (4ld),  Independent 
of  =  —  and  =  —  and  reducing  to  zero  at  w  =  0.  The 

expressions  and  ^e""^®  represent  a  kind  of  motion  between 

two  parallel  walls.  These  solutions  could  he  directly  obtained  by 
in-verting  and  exchanging  the  variables  in  equations  (48)  and 
putting  Wj,  =  nip  =  0. 
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If  n-:3*oo,  eauation  (44)  shows  that  E  ^ Similarly 


the  equations 


1  _  dv^ 


2  ^  1  ^  d(l/m"^ 

“  T_oo^  <i(l/w2) 


coincide  with  l/a^  and  l/p'^  (eauations  (12))  and  can  he  explicitly 
calculated  hy  eq.uation  (l6).  Then  equations  (46)  show  that 


_  -yr  ““1  “"1  __  X. 

^00  =  4.  “  ~  ~  ® 


where  X,  =  log  w  —  /;[*  —  (d  log  m/d  log  w)  d  log  w  coincides 

with  the  value  deduced  hy  integrating  the  dX.  given  hy  equation  (11 ) 

X, 

0 

and  determining  the  constant  of  integration  so  that  =1  at  w  =  0. 

w 

Hence  e^  Is  the  subsonic  asymptotic  value  of  and 

for  n  ^  00. 
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The  solutions  of  equations  (4l)  Just  discussed  are  hound  to  the 
solutions,  equations  (42),  hy  simple  relations.  For  Instance, 


^n(^r) 


Cn(w,m)  +  En(wp)Sn('w,m) 


in(wr) 


”  -  Tn(wr)Sn  -jCT 


as  can  he  yerlfled  hy  controlling  the  identity  of  the  conditions 
at  w  =  Wj,.  Hence  as  w ->  0  and  |Cj^|  and  jSj^|_>co 


11m  Cn(w,m)  _  nm 
v-^0  Sn(w,m)  v-»  0  Cn(m,v)  ” 


C  fii 
lim  ^\m'’v 


sJii 


n  1 

MmV, 


lim  J^yv^m 
n  (II 


The  following  interesting  expansions  are  deduced  hy  applying 
equations  (53)  to  the  solutions  in  closed  form  obtained  for  n  =  +1  : 


=  Wj.C2(w,m)  +  iii^S^(w,m) 


l/w„  "^1 


d(l/w3^)  =-is^(w,m) 


(and  the  corresponding  expansions  with  the  variables  Interchanged  or 
Inverted) . 


1+4 
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FinaJ-ly,  it  should  he  observed  that  differentiating  (for  instance) 
equation  (4la)  with  respect  to  w  yields 


1  ^ 

d 

m  dw 

L  <i(l/m) 

xd-w  / 

=  0 


This  equation  is  of  the  same  general  form  as  equations  (4l)  with 

2  2 

only  one  of  the  variahles  inyerted  and  1  —  n  In  place  of  n  ;  it  can 
therefore  he  treated  in  the  same  vay  as  equations  (4l). 

Hence  two  particular  solutions  similar  to  equations  (42) 


can  he  defined  through  the  integrals  given  hy  the  formula  (26) 

hy  simply  replacing  the  present  variahles*  For  n^  >1  the  second 
series  is  imaginary  and  must  he  divided  hy  i  to  obtain  a  real 
solution.  The  general  solution  for  dX^/dw  is  given  hy  a  linear 

combination  of  these  solutions^  and  the  general  solution  for  is 

obtained  hy  integrating  and  adding  an  approximate  constant.  Now^,  this 
must  coincide  with  the  one  in  terms  of  C^(v^m)  and  S^(w^m)  .  It  is 

then  easily  derived  that 


C 


mp 

m” 


Cj^(m,v) 


m  n 
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and 


Analogous  relations  obtained  by  Intercbanglng  and  Inverting  the 
variables  also  hold.  For  n^  =  1  these  relations  give  eq.uations  (54) 
as  a  particular  case. 


Other  interesting  relations  can  be  deduced  In  the  same  way 


Infinite  series  in  the  exponential  set.—  The  series  in  0^e 
and  have  been  used  (as  series  in  and 


could  be)  by  many  authors  In  the  case  of  the  normal  Isentroplc  law 
(references  5^  6^  9^  sxid  12). 


They  seem  to  have  their  natural  field  of  application  in  the 
problem  of  two-dimensional  gas  Jets^  as  Chaplygin  first  showed  in  his 
classical  memoir. 


The  application  to  flows  around  bodies  seems  to  be  more 
difficult,  especially  for  flows  with  circulation.  The  difficulty 
arises  first  from  the  presence  in  the  holograph  plane  of  a  Singularity 
at  w  =  Wco  and  from  the  ensuing  necessity  of  employing  more  than  one 
series  development  in  the  exponential  set  with  different  sequences 
of  n  (as  appears  already  in  the  incompressible  case)  with  added 
eventual  terms  in  other  sets,  and  of  insuring  that  the  different  series 
are  the  continuations  of  each  other. This  can  be  achieved  (although 
in  a  not  very  simple  way)  by  putting  the  condition  of  continuity  of  the 
solutions  and  of  their  derivatives  on  the  the  transition  curves  (often 
circles),  as  has  been  done  by  Tsien  and  Kuo  (reference  12)  and  as  the 
author  himself  has  done  in  an  unpublished  work  in  a  somewhat  different 
way,  but  it  is  believed  that  the  main  obstacle  to  this  method  arises 
from  the  difficulty  of  insuring  that  the  body  will  have  a  closed 
contour  when  a  circulation  is  present.  In  fact  if  the  so-called 
"natural**  series  (that  is  a  series  having  the  same  coefficients  as  in 
a  chosen  incompressible  case  with  circulation)  is  used  in  one  part  of 


l^This  difficulty  is  avoided  in  the  method  by  Bergman  (reference  6), 
which  uses  a  different  type  of  expansion  and  uses  the  series  in  the 
exponential  set  only  as  eventual  auxiliary  series. 
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the  hodograph  plane  (so  that  a  basic  series  from  which  the  coefficients 
of  the  other  series  -will  he  deduced  hy  the  f oresaid  method  is  obtained) , 
the  resulting  body  •will  be  closed  only  in  the  limiting  Incompressible 
case.  In  a  tentative  method  the  author  has  tried  to  obtain  the  closing- 
up  of  the  contour  by  talcing  the  coefficients  of  the  basic  series  as 
simple  functions  of  E^Cw^)  (or  of  the  other  quantities  in  equations  (43)^ 

Wj,  generally  coinciding  with  Wo()  containing  an  arbitrary  parameter. 

These  simple  functions  reduce  to  the  coefficients  of  the  "natural” 
series  when^  for  vanishing  w^,,  Ej^^  becomes  unity .  The  arbitrary 

parameter  is  then  so  determined  that  the  contour  closes  up.  However, 
because  of  the  necessity  of  using,  to  express  this  condition,  different 
series  connected  by  intricate  relations,  this  method  seems  to  be  very 

complicated.^^  In  the  method  studied  by  the  author,  series  of  the  kind 
given  by  equations  (31),  (32),  and  (33)  (that  can  be  put  in  simple 
relation  with  series  in  the  exponential  set)  coiold  be  used,  especially 
for  the  condition  in  transonic  and  supersonic  regions  where  the  series 
in  the  exponential  set  cease  to  be  useful. 

It  is  worthwhile  to  mention  here  that  the  demonstration  of  the 
convergence  of  the  series  in  the  exponential  set  (that  Chaplygin  first 
deduced  in  the  hypergeometrlc  case  in  a  somewhat  complicated  way)  can  be 
obtained  very  simply  and  under  very  wide  assumptions  for  m(w)  by  using 
the  properties  of  the  functions  defined  by  equations  (43) .  Taking,  for 

instance,  the  series  in  =  (Xne^0)“,  it  is  immediately  seen 

that  Xq  >X^  for,  by  equations  (47)  and  (5l)^ 

« 

^  =  (d  log  m/d  log  =  (1  —  <  1 

^1 


and 


log  - 


[l  -  (E^i)]  d  log  w  >  0 


^^The  problem  of  the  closed  contoixr  has  been  solved  in  a  very 
elegant  way  by  Ll^thlll  (reference  8),  who  discovered  a  very  simple 
development  converging  in  all  the  field  (subsonic  and  transonic)  and 
gave  the  conditions  for  the  closlng-up  of  the  flow  behind  the  body. 
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Moreover,  for  plausible  m(w). 


Rl  =  -£- 

^  PO 


R„=  (1-m2)V2^ 

PO 


dRi  dE 

are  decreasing  functions  of  -w,  so  that  -  <  0  and  — -  <  0.  Taking 

dw  dw  ^ 

nov  a  value  of  n  greater  than  1  it  is  immediately  seen  from 
eq.uations  (kk)  that  if,  for  some  value  of  v,  ^ 

then  >0,  and  if  then  — ~  V  >o.  Hence 

CLW  dw 

if  one  of  the  tvo  conditions  is  verified  for  some  value  of  w,  then  for 

decreasing  w  the  value  of  R^  will  diverge  more  and  more  from  the 

value  of  Rj^  and  of  R^,  so  that  it  cannot  he  equal  to  unity 
at  w  =  0.  Hence,  if  R^  =  1  at  v  =  0,  for  other  values  of  w, 

R^  <  Rn  <  Rp  and  from  equations  (46)  <  Xj. 

The  following  limitations  for  %  =  X^^  are  found  hy  use  of 
equations  (47)  and  (52),  for  n  >  1: 


<Xti< 


Naturally,  these  limitations  hold  only  for  real  values  of  R  and  X 

CO  ^ 

hence  for  M  >1.  In  the  same  way  it  is  proved  that,  for  n  >1 


m^<  ai^<  e^ 


The  analogous  demonstration  for  0^^  and  ijr  requires  the  assumption 

XV  4-  m  2  d(l/w2) 

that  T_ - must  he  an  increasing  fimction  of  w.  This  does 

d(l/in2) 
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not  seem  to  te  too  restrictive  a  condition^  for  T_^  must  in  any  case 

te  unity  at  w  =  0  and  Infinity  at  M  =  1.  According  to  this 
assumption  it  can  "be  directly  verified  that 


mw 

d(l/m2)1 

dw 

/ 

1 + ^  “ ) 

\d  log  w/ 

—  m  —  — 

1  +  ^“8  “f 

\d  log  w/ 

2  dw 

_d(l/w2)J 

>0 


and 


-1 

1  mw^  ^ 

1  iii3 

^  d(l/m2) 

dw  ' 

(m  +  vj 

(m  +  w)^ 

d(l/w2)_ 

By  integrating  these  inequalities  between  0  and  w  it  follows  easily 
that 


the  sign  of  equality  holding  only  at  w  =  0.  Hence  equations  (49) 
and  (51)  show  that,  excluding  w  =  0,  1  <  T_2  <  T  and  the  equation 

corresponding  to  (44)  shows  that  like  T_^  is  an  increasing 

function  of  w.  How  for  -n  <  — 1  from  the  equation  in  by  a 

reasoning  Identical  to  that  developed  for  En,  it  can  be  proved 
that  1  <  T_i  <  T_  <  T_^  and  that  -  >  >  Y_^  >  Y_^.  Hence 

for  the  following  hoimds  hold^  for  n  >  1: 


n 


Wj  dm^^ 


n 


<  e 


nX 
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Similarly,  for 


<  tn  <  I  "^1  '^^ij  <  ^ 


The  convergence  of  the  Chaplygin  (or  other)  series  can  be 
immediately  demonstrated  by  means  of  these  bonds. 

For  negative  Integral  (eventually  half-integral)  values  of  n,  for 
the  reasons  discussed  in  the  section  entitled  "The  exponential  set  "  the 
condition  that  E  and  T  are  unity  at  v  =  0  is  not  sufficient  to 
determine  the  solution.  Hovever,  it  is  possible  to  give  supplementaty 
conditions,  which  are  omitted  for  brevity,  such  that  the  resulting 
solutions  of  eq.uatlons  (4l)  may  be  used  to  construct  series  converging 
in  all  the  subsonic  holograph  field  exterior  to  a  given  circle. 

A  new  set.-  Here,  only  briefly  mentioned,  is  a  different  set  of 
solutions  of  equations  (23)  in  closed  form.  It  has  been  observed  in 
the  section  entitled  "The  exponential  set"  that  the  solutions  0^ 

and  given  by  equations  (5O)  cannot  be  derived  from  the  solutions 

of  the  equations  in  and  for  n^  =  1. 


Conversely,  if  for  given  G  =  0j^e-in0  and  G  =  of  the 

e33)onentlal  set  values  of  F  and  F  are  deduced  by  means  of 
ecluatlon  (25)^  or  the  analogous  eq_uatlon  for  5^^  the  solutions  obtained 
are  still  of  the  corresponding  exponential  set  for  all  values  of  n 
but  1. 

In  this  case,  writing  instead  of  G  =  and  "5  =  t  the 

more  general  formulas 


and 
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corresponding  to  equations  (48),  it  follows  from  equation  (25), 
with  ©p  =  0,  that 


are  new  solutions  of  the  corresponding  equa,tlons  In  (23),  which  are  not 
included  In  any  of  the  sets  already  discussed. 

Then  if  the  varlahles  are  Inverted  and  interchanged  and  the  sign 
of  9  is  changed,  it  is  seen  that 


satisfy  the  corresponding  equations  in  (23);  hence  the  respective  real 
nn(i  imaginary  parts  satisfy  the  equations  for  0  and  i  . 

Wow  applying  again  equation  (25)  to  the  last  expressions,  other 
solutions  of  the  equations  in  F  and  'J'  are  found,  from  which  hy 
changing  again  the  varlahles  new  solutions  G  and  are  found.  Hence 
it  is  seen  that  the  repeated  application  of  the  described  process 
generates  a  new  set  of  solutions  F  and  F,  and  G  and  “fr. 
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The  first  terms  of  the  corresponding  Incompresslhle  set  of 

which  can  he  easily  obtained  hy  repeated  application  of  the  formula  to 
which  eq.uatlon  (25)  reduces  for  m  =  w 


(where  V  =  we 


and  successive  Inversion  of  Y,  are 


Y/2 


THE  GENERALIZED  POTENTIAL  FUNCTIONS 


•  In  the  preceding  sections^  it  has  been  shown  that  the  symmetrical 
form^  obtained  by  making  the  velocity  and  mass  velocity  ^and  not  the 
connection  between  them)  appear  explicitly  in  the  holograph  equations, 
gives  rise  to  an  interesting  general  treatment  of  these  equations.  It 
has  been  seen  that  the  complex  functions  F  and  G^  of  which  the  real 
parts  are  x  05  connected  by  symmetrical  relations  to  the 

functions  T  and  having  cd  and  Y  as  imaginary  parts.  In  this 
section,  it  will  be  shown  that  all  these  complex  functions  can  be 
deduced  by  simple  differentiations  from  a  unique  function  called 
the  generalized  potential  function. 


52 
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A  second  generalized  potential  function  $  is  also  introduced, 
vith  interesting  properties. 

Let  the  complex  velocity  and  mass  velocity  be  defined  by 


and  observe  that  they  are  bo-und  by  the  condition  that  their  ratio 

V  P 

—  =  —  must  t©  real  and  equal  to  a  prescribed  function  of  or 

V  Po 

that  w  must  be  a  prescribed  function  n(w/v)  of  the  relative 
density.  Vhen  these  conditions  are  satisfied^  the  two  moduli  of 
equations  (55)  will  be  connected  by  IWI  =  m(|y|).  Hence  egiiatlons  (55) 
can  be  written 


so  that 


“  -  <}) 


and 
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It  is  now  seen  that  if  W  and  V  are  not  hound  hy  the  foresaid 
conditions  and  are  independent^  these  relations  can  he  used  to  define 
generalized  complex  values  of  ^(YjW)^  and  m(Y_,W)  and  there¬ 

fore  of  all  the  related  q.uantlties.  Hence  the  equation 


1  _  m2  =  m  ^  ^  _  d.  log  P  ^  ^  _  d  log  (W/y)  , 

d  log  V  d  log  w  d  log  n(W/Y) 

defines  a  complex  Mach  nmher  M  =  M(W/y)  whlch^  if  the  physical 
conditions  concerning  Y  and  V  are  satisfied,  reduces  to  the 
real  M(p/pq)  and  can  therefore  he  immediately  deduced  without  the 
help  of  equation  (56)  hy  slirqpiy  replacing  in  the  expression  M(p/pq) 
the  real  varlahle  p/Pq  hy  the  generally  complex  varlahle  WyV. 

Now  the  holograph  equations,  considered  for  complex  values  of  the 
varlahles,  can  he  transformed  hy  taJclng  Y  and  W  as  new  independent 
variables.  From  equations  (55) j 


A 


and 


^  =  e-i® 

Sw 

SY  _  —10  ^ 
Sm  "  ^ 


Sw  dw 


Sm 


e-i® 


de 


5^ 
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so  that 


m 


■N 

=  wf—  + 

Sw  VSv  dw  SW/ 


L.  =  V  ^  > 

^  ^  \dm  Bv  Sw  I 


s  a  3 

1—  =  V--+W  — 

a©  av  aw 


(57) 


Hence  eq.uatlons  (23a)  and  (2313)  136001116 


^(fv  +  ^  ^w) 


+  Hw  ^W)  “  '^T  + 


and 


where  subscripts  denote  partial  derivatives.  These  are  satisfied  If 


or 


"  ®’w 


^  dm  ^ 


=  dj£^ 


(58a) 


(58b) 


d  log  w 


W 


(58c) 
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The  first  of  these  equations  can  he  satisfied  if  a  function  $  is 
introduced  such  that 


F  = 


(59) 


and  the  second  will  also  he  satisfied  if  is  a  solution  of  the 
equation 


=  (1  -  M2)¥2$y^ 


(60) 


p 

where  1  —  M  is  the  function  of  V/V  defined  hy  equation  (56). 

The  meaning  of  equations  (59)  and  (60)  is  the  following.  If 
a  ®(V,V)  satisfying  equation  (60)  is  known,  and  after  Oy  and  are 

calculated,  the  ri^t  values  of  V  and  W  are  introduced  (that  is,  sudi 
that  V/V  is  real  and  IWI  =  iii(lvl))^  then 


^  —  R.P.($y) 


CD  =  I.P.($y) 

will  he  solutions  of  equations  (8) . 

Observe  that  if  in  equation  (60)  the  two  complex  variables  are 
replaced  hy  the  real  variables  9  and  w  (and  m(w))  hy  the  inverse 
relations  of  equations  (57)j  the  resulting  equation  (in  a  complex  $ ) 
remains  unchanged  by  exchanging  w  and  m.  The  functions  F  and  F 
(and  X  and  cd)  can  he  deduced  hy  relations  containing  the  derivatives 
of  t  with  respect  to  w,  m,  and  10. 
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Eq.uations  (23c)  and  (23d)  can  "be  treated  in  the  same  vay  and  the 
result  is  found  that  if  j(l/7,l/W)  is  a  solution  of 


(1  _m2) 


If 


fe  if 


then  the  functions 


will  he,  upon  auhstitution  of  the  ri^t  values  of  ¥  and  V,  solutions 
of  equations  (23c)  and  (23d),  so  that  the  functions 


0  =  R,P.| 


Ilf  =  I.P. 


will  he  solutions  of  equations  (6). 


8 
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Naturally  0  and  xlr  can  also  be  related  to  <I> .  From  equations  (24), 
vlth  the  help  of  equations  (57)  and  (58),  it  follows  that 


\ 


w  w  ''*vw  V  W/ 


G  =  +  —  0^  -  ^1—  +  -jj- 


(63) 


Hence,  comparison  with  equations  (62)  shows  that  to  within  an 
unessential  constant. 


The  physical  coordinates  are  deduced  by  the  followlnfr  relations 
which  may  be  obtained  from  equations  (30)  by  use  of  equations  (55),  (57), 
and  (58): 


S  +  i5f  =  e^^(^<&  +<5  ) 

\W  W  VW/ 


S  +  in  -  e-‘«  ^  « J 


7 


(65) 
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After  assigning  to  W  and  V  their  correct  yalues,  separating  the  real 
from  the  Imaginary  part,  and  using  eq.uations  (2),  it  can  he  shown 
explicitly  that 


^  =  ®vw  ^  ^  ^  ('^WW  ^'^ww) 


which  can  also  he  written 


V  i 


By  differentiating  this  equation,  there  can  he  deduced  an 
expression  of  dz  which  must  coincide  with  the  expression  (2),  that  is 


with 


when  d0  and  dt  are  obtained  hy  differentiating  the  expression 
derived  from  equations  (63) 


2  S 

^  3v\- 


il.P.  W 


It  can  he  shown  that  the  agreement  exists  if  4  is  a  solution  of 
equation  (60),  and  if  W/V  is  real.  It  seems  therefore  that  the 
cLparison  of  the  values  of  z  calculated  directly  from  eiyation  (66), 
or  deduced  hy  integration  of  equation  (67)  (using  equation  (68)),  may 
constitute  an  interesting  check  on  the  accuracy  of  approximate  methods. 
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Oteerve  "that  3>  =  +  c^W  is  a  particular  solution  of 

©q.uatlon  (60),  the  simple  meaning  of  which  is  that  C2_  constitutes  an 
additive  constant  of  m  and  of  -4,  C2  the  same  for  X  and 
and  c^  (complex)  represents  a  general  displacement  of  the  origin  of 
the  physical  coordinates. 

Particular  sets  of  solutions  of  equations  (60)  and  (61)  are  easily 
deduced  from  the  sets  studied  in  the  preceding  sections  ty  integrating 
the  relations  (59)  and  (62)  and  introducing  the  expressions  of  w,  m, 
and  9  as  functions  of  V  and  W, 

In  the  Incompressihle  case,  V  =  W,  1  -  =  1  and  equations  (60) 

and  (61)  reduce  to  identities  satisfied  hy  every  function  of  T.  The 
written  relations  reduce  then  to 


F^(V)  =  <Di‘(v) 


%(v)  =  (V)  =  vei"(v)  -  ^ 

^  dY  \  V  / 

zi(y)  =  Fi*(Y)  =  $i”(Y)  =  e^®(s  +  iU)  =  w 


If  the  solution  corresponding  to  the  incompressihle  flov  around  a 
given  body  is  knovn  in  the  physical  plane^  can  he  deduced  from 

these  relations.  The  profile  of  the  ho(3j  can_,  for  Instance,  he  defined 
hy  a  relation  S  =  P(K)  hetveen  the  suhnorraal  and  the  nonnal  for  'tir  =  0 
(and  hy  giving  this  value  to  the  corresponding  streamline). 

Hence  for  I.P.(Gi)  =  0,  that  is,  for  I.P.(V$^”)  =  it 

follows  that 


wS  =  R.P.(y%")  =  wpji 


6o 
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The  correspondiiag  relation  for  the  compressible  case, 
can  be  easily  deduced  from  the  preceding  formulas,  is 


for  I 


p  w2a 

5w\v  w  / 


=  0 


wMcli 


after  the  correct  values  for  V  and  W  are  introduced.  The  function  P 
will  be  the  same  in  both  cases  if  the  profile  is  imchanged. 

The  eq.uati’ons  (60)  and  (61)  are  of  the  same  general  kind.  It  is 
possible  to  pass  from  the  one  to  the  other  not  only  by  a  substitution 
like  eq.uation  (64)  but  also  by  siii^jly  putting 


$  =  YT#* 

>  (69) 


J 


since,  as  can  be  immediately  verified,  must  satisfy  equation  (61) 
and  f*,  equation  (60).  The  equations  may  be  transformed  in  many  ways 
by  changing  the  two  Independent  variables.  One  of  these  transformations 
is  obtained  by  taking  as  new  independent  variables  (K  being  a  constant): 


e  1  1 

_  1  1 

^  V  VyfK 

J 


(TO) 
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The  equation  (6l)  in  $  (or  41*)  is  thus  transformed  into 


In 


(71) 


where  the  first  factor  of  the  right-hand  side  is  a  function  of  W/Vj 
hence  of  y\/\,  since 


]L  =  1. 1  +  (n/l) 

V  VE  1  -  (tj/I) 


A  similar  transformation  may  h©  performed  on  eqiiation  (60)  hy 
putting  I  =  V  +  W  and  =  V  — 


Other  interesting  transformations  are  obtained  by  talking 
I  =  log  V  +  €{V/V)  and  tj  =  5(w/v)  and  choosing  in  different  vays 
the  functions  e  and  5.  In  this  case  the  variability  of  the 


T|  variable  can  be  restricted  to  the 
obtained  by  taking^  for  instance^  € 


5  =  1 

case 

same 


and  5  =  w(V/v) 


taking 


^ since 

is  the  one  for  vhlch  I  =  X,  —  10 
as  in  equation  (ll)  or  equation 

de  -1 


real  field.  Particular  cases  are 
=  0,  e  =  log  (V^), 

An  Important 


log  V 

d  log  (V/W) 

and  T[  -  Xy  where  X  is  the 

(52);  this  is  obtained  by 

d6  -M^ 

and  -  =  - 


If 


d  log  (y/W)  (1  +  VI  -  M^)  d  log  (y/W) 

this  transformation  is  performed^  an  equation  in  0*  and  f  (or 
in  ^  and  $*)  is  deduced  that  may  be  used  to  obtain  directly 
solutions  of  the  kind  obtained  by  Bergman  (reference  6)  and 
Lighthill  (reference  8) . 


one 


Finally^  let  equation  (60) 
law.  In  this  case 


be  written  with  the  actual  isentroplc 
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Hence^  consideration  of  the  observation  following  equation  (56)  and 
substitution  of  the  resulting  value  of  MCv/v)  in  equation  (60)^ 
yields: 


WW 


There  is  a  corresponding  equation  for  equation  (61).  For  real  W/V  the 

P 

corresponding  values  of  1  —  are  shown  in  figure  5  ^or  some  values 
W  P 

of  7  with  —  =  —  as  abscissa.  It  is  interesting  to  observe  that 
V  Pq 

if  7  is  in  the  actual  range  for  gases ^  its  value  does  not  seem  to 
affect  to  a  great  extent  the  shape  of  the  curves,,  especially  in  the 

subsonic  range.  The  factor  1  —  is  linear  in  V/V  for  7  =  0  and 
in  V /W  for  7  =  2;  for  other  values  of  7  in  the  actual  range  it  is 
not  far  from  a  straight  line  in  the  subsonic  range. 


The  curve  7  =  — 1^  that  is 


p 

1  -  -  — 


is  also  represented  in 


the  figure^  and  corresponds  to  Chaplygin's  approximation.  The  Karman- 

[2  =  ssd 

,.2 


Tsien  approximation  corresponds  to  1 


where  the  constant  K 


(see  equations  (1T)  and  (21))  Is  so  chosen  that  at  Infinity  (p  =  p,^) 
1  -  will  take  the  value  1  -  given  hy  the  true  law. 


THE  CHAPLYGIN-K^M^-ffSIEN  CASE 


For  the  Chaplygin-Kirman— Tsien  approximation  the  right-hand  side 
of  equation  (7I)  (or  equation  (69))  becomes  zero  so  that  the  general 
solutions  of  equations  (60)  and  (61)  are 


f  =  e2(  I)  +  e2(Ti) 


(72) 
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■with  ©2,  and  f2  arhltraxy  f-unctlona  of  the  variahles  | 

and  T)  defined  hy  eq,uatlons  (70) .  The  value  K  =  1  corresponds  to 

Pq^  /  / 

Chaplygin'' s  approximation;  K  =  (l  -  — p,  to  the  Karman-Tslen 

Poo 

approximation.^^  As  It  has  "been  observed  In  the  section  entitled 
"Approximate  Methods"  it  can  "be  convenient  to  chose  K— values  'between 
the  two. 

The  solutions  (72)  can  he  interpreted  in  two  different  ways. 

First  they  can  he  regarded  as  the  exact  solutions  of  the  corresponding 
eq.uations  for  a  gas  satisfying  the  ideal  law  (18) . 

With  the  use  of  eq.uatlons  (7O)  for  |  and  T|,  this  law  can  he 
■written  lilhl  =  -E;  or  for  real  W/Vj  hence  for  real  |  /r\. 


IIt)!  =iT)  =  'nT  =  -fl  (73) 

This  relation  allows  the  expression  of  f  (and  as  the  sum  of  two 

arhltrary  functions  of  |  and  f  or  of  t)  and  t). 

In  the  second  Interpretation,  equations  (72)  are  considered  as 
approximate  solutions  of  the  equations  for  the  actual  law  of  gases.  In 
this  case 


U  T,|  =  iTf  =  TiT  =  i  (1  -  ^ 


P% 


(74) 


is  no  longer  a  constant,  hut  is  a  function  of  w. 


^^As  observed  in  footnote  4  the  actual  presentation  of  this  method 
is  more  coherent  than  the  usual  one,  as  the  constants  H  and  K  of 
formula  (18)  are  here  deduced  for  a  single  ref erence . condition;  namely, 
the  infinite  point. 


61^ 
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Naturally  the  last  expression,  and  the  relation  (73) j  not 
restrict  the  Independency  of  |  and  t|;  as  they  must  he  used  only 
after  all  the  foimial  deductions  from  equations  (72)  of  the  following 
kind  have  been  performed,  l^om  equations  (72)  and  the  application  of 
equations  (59)^ 


F  =  $  =  Y 

¥ 


f  — h=-  '  +  fp  +  — ^  fp') 
\  ^  Wn/k  ^  ^  wVf  ^ ' 


F  =  $, 


and 


Y 


W 


2fi  -  +  2f2  -  t)f2’ 


According  to  eq,uation  (64),  the  last  expression  coincides  vith 
hence,  hy  equations  , 


e-,  = 


®2  “  ^^2  '^^2' 

From  equations  (62), 

G  =  0  +  1^  =  If 3^"  -  +  Tif2"  - 

and 

S'  =  '^  +  It  =  i  (if 3^"  _  f3^»  -  Tjf2"  +  fg') 
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(77) 
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This  relation  gives  the  physical  coordinates  as  a  function  of  the 

.  fl  1  1  19 

holograph  coordinates  w  and  9  vhen  I  -  -  + 


and  Tj 


law  m(w)). 


.(»)n/K 


introduced  into  it  (with  the  assigned 


Differentiating  equation  (77)  yields 


2  dz  =  |2fi"‘dl  +  +  |lTj|(fi"'dl  +  f2*"dTi)  +  (fl'*  +  fa") 


But  dz  is  also  given  by  equation  (2),  which  may  be  transformed  in  the 
following  way 


=  I  d(0  +  1  +  r\  d(0  -  1  \f&) 


or  (see  equation  (75)) 


2  dz  =  |2fi'”d|  +  Tj^f2'**d.Ti  +  |N|(fi”'‘^^  +  f2"’d.n) 


Hence,  in  accordance  with  the  observation  following  equation  (68),  these 
two  expressions  for  dz  coincide  only  if  llnl  is  constat;  that  is, 
if  the  law  connecting  I  and  tj  is  the  law  (73)  ,))^  -i+^+hR 

same  law  as  the  one  for  which  the  equations  (60)  (61) 

solutions  (72).  This  is  what  is  done  in  the  Karman-Tslen  method. _  If, 
on  the  contrary,  the  alternative  interpretation  of  equation  (72)  is 
adopted,  the  coincidence  ceases  to  exist.  In  this  case,  if  the  la 
connecting  I  with  is  the  exact  gas  law,  the  error  term  between  the 

two  values  of  ■fclia't  Is^ 


+  f^")  dllTjl 


with  IIbI  given  by  equation  (74),  may  be  reg^ded  as  a  measure  of  the 
approxlma'tion  of  "blie  approximate  solution  (72)  • 
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Observe  that  If  the  lav  (73)  is  taken,  hence  if  if  =-^, 

equation  (75)  shows  that  the  general  solution  for  0  +  i'®f  is  a 
fimction  of  I  alone  (containing  the  constant  H),  that  is,  it  is  a 
function  of  X  —  10,  where  X  is  given  hy  equation  (19),  a  result  in 
accordance  with  the  Chaplygin's  monogenelty  conditions. 

The  expression  (77)  for  z  can  he  also  written  hy  using  the 
functions  g^^  and  g^  instead  of  f^  and  f^  so  that  (see 

equation  (76)), 


2z 


|2fi"'dl  + 


hl|(fi"  +  f2”) 


(78) 


If  the  law  (73)  ia  assumed. 


this  equation  reduces  to 


2z 


=  J  +  dgg)  +  J +  dgg) 


(79) 


where  g(|)  =  g3^(|)  +  g2(-^[/|)  is  the 
this  case.  In  the  incompressible  case 
2 

with  |-  =  V .  Hence  Tsien’s  formula  z 
immediately  obtained. 


complex  potential  0  +  l>fE^  in 

H  =  0^  K  =  and  2z^  =  ^  dp; 

\2  ^ 


=  Zn-  — 


dz.-  is 
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It  Is  well-knowi  that  in  the  case  of  the  flow  round  a  body  Tsien's 
formula  generates  closed  profiles  only  if  circulation  is  absent.  Many 
authors  have  studied  eztentlons  of  the  method  to  the  case  with 
clrcxilation.  BerSj  Germain^  and  Leray  (references  13  and  l4)  have 
followed  a  first  way;  Lin,  Germain,  and  Gelbart  (references  15,  1^, 
and  7)  a  second  way;  here  a  third  way  of  constructing  flows  aromd 
closed  profiles  with  circulation  will  be  shown,  based  on  the  subdivision 
of  g(|)  into  gi(|)  and  g2(l)*  Let  gi(‘V„/7(zi) )  =  gi(l/l(»)  "be  the 
complex  potential  of  an  Incompressible  flow  around  a  closed  profile  with 
a  circulation  F.  Then,  since  for  |  near 

(l/IJ^  =  1  +  k[(l/lco)  -3+0  |[(|/U  -  ij^j 


it  follows  that 


5  dgi  =  2  ^dzi  =  0 


/  dgi-s/idzi— ^ 


=r 


where  the  integrations  are  performed  along  any  contour  in  the  physical 
plane  enclosing  the  profile,  or  around  the  corresponding  contour  in 
the  I— plane  enclosing  (for  simplicity,  suppose  and  real). 
Wow  generally. 


l^dzj^  =  k| 


-/(I- 


.  X  ,  Ic  .  k+1  „ 

ljdzi=--|^  r 


SO  that 


=  (8.) 

ajid  the  value  of  z  given  hy  equation  (79)  for  a  compressible  flow 
with  g  =  gj  is  not  one-valued.  But  if  g  =  g^  +  g2 

given  gj^^  a  g^  can  be  determined  in  such  a  way  that  the  corresponding 
residual  terms  in  equation  (79)  will  compensate  the  value  of  the 
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last  vrltten  Integral.  If  the  value  of  the  circulation  must  remain 
unchanged,  gg  must  he  one-valued.  This  condition  is  obtained  very 

simply  hy  taking,  for  instance. 


The  constant  h  can  now  he  determined  so  as  to  obtain  the  said 
compensation.  From  the  identity 


and  from  equation  (8l)  It  follows  tliat 


r 

n  +  r 


21i? 


n-1 

00 


-^r 


Hence  for  real 


loo  and  H  =  -I„ti„o 


<162 


82  =  "^(^00  -  Or 


If  eq.uation  (79)  m-ust  he  one-valued,  the  last  quantity  must  he 
equal  to  (see  equation  (82)) 
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Hence  the  equations 


2TJoo 

h  =  - 

too  -  nco 


and 


gg  = 


(83) 


will  satisftr  the  foresaid  .condition  of  compensation. 

Ohserve  now  that  the  constant  factor  does  not  depend  on  n.  Hence 
the  condition  is  also  satisfied  if  is  replaced  hy  P(  loo/t)/P(l), 

where  P  is  a  polymonial  in  or  an  infinite  series  converging  in 

all  the  domain  of  variation  of  Hence  the  fnnction 


0  +  =  g(l)  =  gi(l) 


too  ^^00 


p(Wt)  n 
p(i)  j 


(84) 


where  P  is  an  arbitrary  function  of  loo/t j  Is  analytic  in  all  its 
domain  of  variation,  and  will  generate  a  flow  around  a  closed  profile 
with  circulation  P , 

Observe  that  by  using  eq.uatlons  (70)  for  real  and  the 

expression  for  the  constant  h  is 

h  =  2!1,«s —  =  £22  \/K  -  1 

too  "  V  ^0 
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For  the  Chaplygin  or  E^rman-Taien  values  of  K  there  is  obtained, 
respectively,  the  very  simple  expressions,  both  vanishing  in  the 
incompressible  case: 


h  = 


1 


and 

h  =  \/l  -  M  ^  -  1 
^  00 


Particularly  simple  forma  of  equation  (8lj-)  are  obtained  by 
putting  n  —  0  or  n  =  1  in  eq^uation  (83).  For  n  =  1  an  expression 
is  obtained  which  coincides  with  what  becomes  the  Llghthill  solution 
for  7  =  -1  (reference  8) . 

The  solution  (84)  satisfies  the  condition  of  generating  solutions 
around  closed  profiles  with  circulation  when  |t)  ||  =  Constant,  and 
z  is  given  by  equation  (79)*  Now  if  jt)  ||  is  variable,  as  given  by 
equation  (74),  and  the  equations  (72)  are  considered  as  approximate 
solutions  of  the  exact  equations  (60)  and  (6I),  it  is  still  possible  to 
find  solutions  for  which  z,  given  by  equation  (78),  is  one-valued.  If 
again  gj|^(|/|^)  is  the  complex  potential  for  the  Incompressible  case, 

satisfying  equations  (80),  with  in  place  of  P,  and  if  ^(irf/Tlco) 

is  a  function  for  vhich  the  expressions 
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(analogous  to  equations  (80))  hold  where  the  Integrations  are  performed 
in  the  tJ"  plane  along  the  contour  corresponding  to  the  one  of 
equations  (80),  then  equation  (82)  and  the  analogous  relation  for  g^ 

give 


80  that  the  sum 


will  he  zero  for 


r„  (86) 

^  boo 

Hence  equation  (78)  will  he  one— valued  if  gg  satisfies 
equation  (85)  with  r2  given  hy  equation  (86),  and 


0  +  =  g^d)  +  ggCn) 
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will  represent,  when  the  expressions  (70)  for  I  fl-nfi  t|  are  replaced 
hy  the  exact  connection  "between  m  and  w,  an  approximate  solution  of 
the  exact  equations  with  circulation 


z  being  given  "by  equation  (78),  and  the  error  term 


representing  a  measure  of  the  approximation  obtained. 

Clearly  the  simplest  way  of  satisfying  equations  (85)  and  (86)  is 

that  of  taking^3 


ggClj)  =_^  gj^(T)) 


Natnreilly  this  solution  still  holds  when  IN  =  Constant. 

An  approximate  solution  for  the  transonic  case  Is  now  noted, 
corresponding  to  the  subsonic  Karmiii>~-Tslen  approximation.  If  =  1, 

then  the  Earni^n—Tsien  value  (equation  (21))  for  K  is  zero.  The 


^3if  the  Karmanr-Tslen  value  (equation  (21)) 
which  makes  the  right-hand  side  of  equation  (7I) 


is  adopted  for  K, 
zero  at  infinity,  then 


4 


f 

^  00 

T 


1-  \/l- 

1  +  \/l  - 


and  1  -1? 
I 


2 _ 

1  +  Vl  - 


i 
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corresponding  curve  of  1  —  in  figure  5  reduces  to  the  horizontal 

axis  1  —  =  0.  In  this  case  the  right-hand  sides  of  equations  (6o) 

and  (6l)  are  zero,  and  the  respective  solutions  are 

^  =  fi(W)  +  VfgCW) 

and 

i  =  i  ei(¥)  +  e^Cw) 


where  the  arhitrarj  fmctions  inYolred  are  bound  by  the  relations 


— e-i  =  f*- 


and 


Eq^uations  (66)  and  (68)  then  give 


z  =  iR.P.CWfg')  +  +  I.P.(f2)] 


0  +  iijf  =  -R.P.(f2^’) 


+ 


+  W^fg" 


+  Wfg'  -  f 


The  approximate  gas  law  corresponding  to  1  —  M  =  0  is^  by 
equation  (56)^  m  =  Constant,  therefore  Pq/p  is  proportifpal  to  w. 

Actually  m  has  a  maximum  at  M  =  1,  and  the  distance  between  the 
streamlines  reaches  there  a  minimum.  If  the  approximate  law  m  =  Constant 
were  adopted,  the  distance  between  the  streamlines  would  be  unchanged 
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throu^out  the  field  of  motion,  and  this  would  give  rise  to  difficulties. 
But  If  the  solution  Is  considered  as  an  approximate  solution  of  the 
exact  eq^uatlons  In  the  transonic  field  (even  for  1  but  near  1), 

then  the  exact  law  for  m(w)  can  be  Introduced  In  the  solutions.  This 
approximate  transonic  theory  seems  worthy  of  development. 


September  27^  19^9 
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Figure  2.-  Comparison  of  the  Kam^-Tsien  and  Chaplygin  approximations 

with  the  true  isentropic  law. 
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Figure  3 .  -  Vaxiables  a  and  P  as 


functions  of  X  for  isen tropic  gas . 
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Figure  4.-  VarialDles  and  as  Figure  5--  Comparison  of  the  Chaplygin  and 

functions  of  for  7  =  1.4.  Karm^-Tsien  approximations  of 

1  -  as  a  function  of  p/Pq  vith  the 
exact  relation  for  several  values  of  7. 


